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Abstract 

A theory of dissipative generalized continuum mechanics is presented in the 
framework of weakly nonlocal non-equilibrium thermodynamics. The evolu- 
tion equation of microdcformation is obtained by thermodynamic principles. 
Conditions of generic stability, the linear asymptotic stability of homogeneous 
equilibrium, is derived in a simple but representative case. 



1. Introduction 

There is a kind of stability that is among our basic and minimal physical 
expectations in continua. Explicitely or implicitely we think, that small per- 
turbations of a homogeneous thermodynamic equilibrium state decay. Small 
excitations of arbitrary form are damped and this damping is expected as a 
result of material properties, otherwise the continuum, as a material, cannot 
exist, it destroys itself. Theories should reflect this property, or they are wrong 
models of a natural system. A relatively simple necessary condition of the above 
mentioned expectation for a specific continuum model is generic stability: the 
linear asymptotic stability of homogeneous solutions of the evolution equations. 
If the conditions of generic stability do not require anything beyond thermody- 
namic reasons, i.e. by the concavity of the entropy and by the nonnegativity of 
the entropy production, then one can be sure that the parameters of the model 
family represent material properties, the mentioned basic stability expectation 
is fulfilled. Generic stability may be a benchmark on the viability of theories of 
dissipation. The more the conditions of generic stability are connected to ther- 
modynamic principles, to concavity of entropy density and increasing entropy, 
the higher can be our hope that our dissipative theory is a viable model of a 
natural system. 

The concept of generic stability was introduced for special relativistic dissi- 
pative fluids by Hiscock and Lindblom [l[ . There generic stability distinguishes 
between thermodynamically consistent looking theories with increasing entropy 
production. There are violently unstable ones conditionally stable ones [2] 
and there are those where generic stability does not require anything beyond 
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basic thermodynamics requirements [3| , and in this sense they are similar to the 
nonrelativistic, dissipative Navier- Stokes fluids Q. 

Thermodynamics and stability of some simple equilibria are considered 
closely related also in solids. Sometimes the expectations are high, failures 
are interpreted as fatal and any relation is denied On the other hand, new 
researches about the fundamental connection of the second law and stability of 
thermodynamic equilibrium Q change the failures to challenges that may lead 
to new approaches to old problems [7|. 

In generalized mechanics dissipation looks like a secondary aspect. Mindlin 
introduced a variational principle Q , Eringen and Suhubi used a balance of mi- 
cromomentum with averaging [9( and Germain suggested the principle of virtual 
power [Icj to obtain new evolution equations. None of them was interested in 
dissipation. The situation did not change much since then, the search for the 
proper interpretation and nondissipative material parameters suppressed the 
need of a dissipative generalization [Tlj . However, recent research revealed a 
fundamental connection of thermodynamics and generalized mechanics. From 
a thermodynamic point of view the microdeformation can be considered as an 
internal variable and the evolution equations of Mindlin can be obtained by 
introducing dual tensorial internal variables [l2[ • 



In this paper we obtain the Mindlin-Eringen-Suhubi theory of microdefor- 
mation by weakly nonlocal non-equilibrium thermodynamics. Hence we arrive 
at a dissipative theory from the very beginning. Then we investigate the generic 
stability of the obtained dissipative theory in a reasonably simple, but repre- 
sentative specific case. It turns out that pure thermodynamics conditions are 
not enough, instabilities may arise. 

2. Mindlin theory of microdeformation 

The fundamental kinematic quantities of generalized mechanics are the dis- 
placement u l — x l — X 1 and the micro- displacement u n = x n — X n , where 
x l and x n are spatial X 1 and X' 1 are material position vectors. Then one 
introduces the strain as the symmetric part of the displacement gradient 
£ l i = gU u v — 1(9V + d^u" 1 ), and microdeformation as the gradient of the 
microdisplacement ip^ = d^v! % . With displacements and Cartesian coordinates 
the treatment is frame and coordinate dependent from the beginning. However, 
the indices may be considered as abstract, that do not refer to any particular 
system of coordinates, but conveniently denote the tensorial components [l3j . In 
the following we will use the indices in this sense. We do not distinguish between 
vectors and covectors. Einstein convention of repeated indices is applied. 

Secondary kinematic quantities are the relative deformation ^ — d % u 3 —ip" 1 - 1 , 
the microdeformation gradient (i i: > k — d k ip l: > and the velocity v l . The velocity 
and strain fields are related by 
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The fundamental balance of momentum is: 



gv l - dH lJ = 0, 



(2) 
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where g is the density, the dot denotes the substantial time derivative, the 
symmetric t ij is the stress and d % is the spatial gradient. The balance of internal 
energy follows as usual: 

ge + SY = t ij d j v i (3) 

where e is the specific internal energy and q l is the current density of the internal 
energy, the heat flux. The evolution of the microdeformation is given as: 

S k ^ k -d k n ijk =a ij . (4) 

Here d lk is the tensor of microinertia, ^ k is the double stress and cr lJ is the 
relative stress. In the ideal, nondissipative theory of Mindlin this equation is 
derived from a variational principle and the three stresses are derived from a 
free energy. Generalized mechanics is a first order weakly nonlocal, gradient 
theory in the microdeformation therefore the free energy is a function of the 
strain, the relative deformation and the gradient of the microdeformation. For 
our purposes it is convenient to consider the free energy as a function of the 
microdeformation instead of the relative one. The two representations are not 
equivalent, because relative deformation depends also on the antisymmetric part 
of the displacement gradient. That seems to contradict to the usual formula- 
tions of material frame indifference therefore in the following we consider 
specific free energy as function of microdeformation instead of the relative one: 
w(e lJ \ i[>' lJ , partial k tp l: ')- Then the equilibrium stresses are expressed by the 
following derivatives of the free energy: 

ij _ dw_ ij _ dw ijk _ dw 

&«' w ^ E ~ d(d k ^y [ > 

In case of isotropic materials a quadratic free energy has the following form: 

w{e ij ,i; ij ,d k ip ij ) = 

*( e «)2 + ^iji + |(V>") 2 + | W j + | W + 
Ji^+jj^+f)^ + 

ai d k i>"d j i> kj + a-2,d k i> il d 3 i> 3k + ^d k ip"d k if) jj + a 4 d j ^ ''d k ip ik + 
a 5 d j 'ip ij d k il) ki + ^-d i ^ ij d k tp jk + ^-d k ip ij d k ilJ ij + — d*^ d i ip jk + 
'-i;) 1 ,/■'<■)' i . ' L + ^-d k ip ij d k ip ji + ^-d k ^ 3 d l ip kj , (6) 

where A and \i are the Lame coefficients, b\, &2, ^3,5i,52: &i> ■ On are material 
parameters. Then, substituting the constitutive functions ([5]) into the evolution 
equations @, © and ^ we obtain a closed system. 

It is generally believed, that Mindlin's version of the microdeformation the- 
ory is nondissipative. Actually, without calculating the entropy production one 
cannot say anything about dissipation. The approach of Eringen and Suhubi 
is superior from this point of view, they calculate entropy production. On the 
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other hand their dissipation - by artificially introducing dissipative parts of the 
stress, relative stress and double stress - is questionable. In classical continuum 
mechanics the entropy inequality is a conditional restriction, it is valid with the 
condition of the validity of evolution equations, the classical balances. In case of 
a second order evolution equation like that condition is not compatible with 
energy dissipation formula of Eringen and Suhubi, eq.(5.13) in Therefore in 
the following we suggest a simple alternative approach of the Mindlin-Eringen- 
Suhubi (MES) theory. 



3. Thermomechanics with dual internal variables 

Based on (l5| one may propose a pure thermodynamic approach where the 
microdeformation ip 13 is an internal state variable and introduce also a second 
tensorial internal variable denoted by /3 13 [IH ■ Then we assume, that the entropy 
density is a function of the internal energy e, the strain e %3 , the internal variables 
ip %3 \ j3 13 and the gradient of the microdeformation d k "p %3 . Moreover, we do not 
require a variational principle or a balance of micromomentum, the evolution 
equations of the internal variables ip i3 and (3 lJ are 

¥ j =f\ ^=g lJ . (7) 

Here the right hand side of the equations are to be determined by the second 
law. We expect to obtain the MES theory as a nondissipative special case. Our 
only non straightforward assumption is the following, non classical form of the 
entropy current density: 

Jl = t ~ e d{d k ^) r ■ (8) 

Here T is the temperature, defined by d e s = 1/T. 
The entropy balance can be calculated as follows: 

< gs{e, e' 3 ,^ 3 ,d k ^ 3 ,(3 13 )) + d t J l = 

gd e e + gd e i]S£ lj + gd^jSip 3 + gdgk^ij s(d k tp 13 ) + dpas^ + 
,sq l - gd d i^ksip 3k ^ = 

- (r ■■<)■<■■ - ,vy) + Q d sir,- + ,><),,. sf 3 - <,<)'■ (d d ^ s ) ./"•' - 
Q d a ^sd l ^d k v l + on , ,<r' + d l (±q^j = 

- (*« + e T{d e , lS - d d ^ klS d 3 \p kl )) av + e {d^ts- d k {d aHijS )) p 3 + 

Qdpn sg* 3 + (q* - gd di ^ sf 3k ) 9, i (9) 

Here we have applied the kinematic condition ([T]) and the evolution equa- 
tions of the momentum (J5|), internal energy @ and the internal variables ([7]) 
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as constraints. In small strain approximation the density can be considered 
constant. Remarkable are the nonclassical extensions of the mechanical and 
thermal contributions of the entropy production. Now we change the thermo- 
dynamic potentials and introduce the more traditional Helmholtz free energy 
representation. Then the transformations of partial derivatives are 



ds 
ds 



1 dw 



T dev 



ds 



dif)^ 



1 dw 



dpi 



1 dw 



T dfri 



dd k Tp^ 



1 



dw 



Td{d k ^) 



(10) 



Substituting these formulas into (|9]) we obtain the following temperature 
multiplied form of the entropy production, that we call dissipation: 

< TP S = (t ij - g(d e aw - d d ^ k iwd j ij kl )) <9V - 

g (d^w - d k (d dk ^w)) f 11 ~ •»>. wg ij + (q l + od„ ■ wf jk ) Td^l) 

This is a quadratic form of terms that are products of the constitutive func- 
tions t l i , f l i , g % i and given functions of the constitutive state space. The last 
two constitutive functions determine the evolution equations. In this case we 
can apply a linear approximation of the constitutive functions in order to solve 
the inequality. It is legitimate to introduce thermodynamic fluxes and forces in 
the sense of non-equilibrium thermodynamics: 





Thermal 


Mechanical 


Internal 1 


Internal 2 


Fluxes 


<f = q l + 

Qf jk dgi^ihW 


ty = t« - Q8&W+ 

gd g i^kiwd : >ip kl 


Qf J 


Q9 V 


Forces 


<9 l lnT 




X*J = 
d k d a k. l f, i] w 


Y*J = -d i ,aw+ 
-dps w 



Table 1. Thermodynamic fluxes and forces of generalized mechanics 

Here t l J denotes the viscous stress. In case of isotropic materials the linear 
relations between the thermodynamic forces and fluxes are the following: 



<f = XdilnT 

i« = L\\ kl d {l v k) +L% kl X kl + L% kl Y kl 

efo = Li jkl dV V V +Lg kl X kl + Lg kl Y kl 

ggV = Lg kl d {l v k) + Lg kl X kl + Lg kl Y kl , 

where for all I, J — 1,2,3 the conductivity coefficents Ljj are 

L)f = suS^SM + auS^SM + l u S^5 kl , 



(12) 
(13) 
(14) 
(15) 

(16) 



where the braces ( ) denote the traceless symmetric part of the correspond- 
ing tensor in the related indices S'&SM = (6 ik 6^ 1 + 5 il S jk )/2 - S ij S kl /3 and 
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the rectangular parenthesis [ ] denotes the antisymmetric part as S^S^ 1 ' — 
(# lfe <P' z — 8 ll 8^ k )/2. This kind of decomposition is instructive because symmet- 
ric, antisymmetric and spherical second order tensors are mutually orthogonal 
in the "double dot" product of second order tensors, when taking the trace of 
their product. Therefore the constitutive equations (fll?T) - (fT5)) can be decomposed 
into three parts: the five component traceless symmetric, three component an- 
tisymmetric and one component spherical parts are independent. Moreover, the 
vectorial and tensorial thermodynamic interactions are independent, too. 

The evolution equations of generalized mechanics are obtained when the 
internal variable (3 lJ is formally a generalized momentum of the internal variable 
tp lJ . One can see that applying two conditions. 

1. Quadratic free energy. The free energy of Mindlin ([6]) must be supple- 
mented by quadratic isotropic f3 13 dependent terms: 

w{e l \ ^,d k ^,^) = w M (e ij ,^ ij ,d k ip i: >) + ^(3 ij B i]kl (3 kl . (17) 

and B ijkl = bsS^SM + b a S^ k 8 jl ^ + biS lj 5 kl . 

2. Ideality. This condition switches off some dissipative couplings. For ex- 
ample for the conductivity coefficients we prescribe that L^ kl = L l ^ kl = 
L33 = % i kl and require antisymmetric coupling of the evolution equa- 
tions of the two internal variables by L 1 ^ 1 = —L^ 1 — L^ kl . 

Then taking the time derivative of (fT4)) and substituting $^ from (fT5|) results 
in the following second order evolution equation of the microdeformation ip %0 : 

gpjkl^kl _ Qk^jk = a £ + pjkl L klmn im n + jijkl jUtku f& mn + gk^mnkj ^ 

Here I ijM = L ijmn B mno P L°P kl is the microdensity tensor. 

The evolution equations (U), ([2]), © and (fl8|). together with the dynamic 
constitutive relations (Tr2 |) -([T5 |) and static constitutive relations ((6J, ([3]) and 
(fTT]) form a dissipative generalization of the MES theory except the kinematic 
interpretation of the internal variables. The evolution equations of the original 
nondissipative Mindlin-Eringen-Suhubi theory is recovered if A = and L^ k 1 = 
for J, J = 1,2,3 except L^ kl = ~L% kl = L^ kl ^ 0. 



4. Generic stability in one space dimension 

In this section we investigate the conditions of generic stability in case of 
a single space dimension, when the fields depend on the time r, and a single 
spatial coordinate, denoted by x. The extension to three space dimensions is 
straightforward . 

The one dimensional forms of the relation of the strain and the velocity (fTJ), 
the evolution equations of momentum ([2]) , internal energy ([3]) are the following 

e-d x v = 0. (19) 
gv-d x t = 0, (20) 
ge + d x q = td x v, (21) 



G 



Here Cartesian coordinates are introduced and the first, x component of 
vectors and tensors is used in the reduction: v(t, x) = v x (t, x), q(r, x) — q x (r, x) 
are the x components of the velocity and heat current fields, t(r, x) — t xx (r, x), 
e(t,x) = e xx {t, x) are the 11 components of the stress and the strain fields. 
d x is the derivative by x and the dot denotes sustantial time derivative, For 
example e = d T e + vd x e is the substantial time derivative of the specific internal 
energy and d T is the partial time derivative. The evolution equations ([7]) of the 
internal variables ip(r,x) — ip xx (r,x) and /3(r, x) — /3 xx (t,x) are reduced to 

1> = f, (22) 
$ = g. (23) 

Here / = f xx and g — g xx . The quadratic, isotropic free energy (fTT)) in one 
space dimension is 

w(e, 4, d x *p, (3) = | e 2 + \V + c # + |/3 2 + \ [d^f . (24) 

The convexity of the free energy requires, that the E, b, B, a coefficients are 
positive and Eb — c? > 0. The dissipation (fTTj) is obtained as 

< -d x (lnT)(q + gfdg^w) + (t- g{d s w - do^wd x w))d x v - 

(d^w - d x (dg m ^w)) gf - (dpw) gg. (25) 

The consequent linear conductivity relations are 

9 =q + Qfdd a ipw = -Xd x lnT (26) 

t v = t - (Ee + ap) + a{d x i>) 2 =l 11 d x v + l 12 gX + l 13 gY (27) 

i> =f =hid x v + l 22Q X + l 23 gY (28) 

(3 =g = hAv + l 32 gX + l 33 gY. (29) 

where gX — —gd^w + gd^da^w) = —trip — ce + ad xx ip, and gY = —gdpw = 
—B(3 and E = gE, b = gb, c = gc, B = gB and a = ga. The inequality of 
second law (|25|) requires, that A > and the symmetric part of the matrix lu, 
I, J = 1,2,3 is positive definite. Here the known transport coefficients in Ijj 
are In, the coefficient of viscosity, and L = (l 23 — l 32 )/2, the antisymmetric 
coupling term of the dual variables, from the microdensity I = L 2 B. 

4-1. Homogeneous equilibrium 

We consider the kinematic condition (fT9")l . the evolution equations ([2H)l - 
(|23|) and the constitutive functions in ([26]) - ([29]) . defining q, t, / and g. It is 
convenient to introduce the stress t as an independent variable in addition to 
e,v,s,if) and /3. Let us denote the homogeneous solution of these equations by 
eo, vq, £q, to, i/jq, Pq. The space derivative of these fields vanishes, they are func- 
tions only of time. Therefore from ([TJ follows Eq — const, and we choose the 



7 



homogeneous equilibrium strain zero. From the momentum balance we ob- 
tain, that wo = const., and from the energy balance ([3]) follows that eo = const.. 
The homogeneous equilibrium values of the internal variables are determined 
by the differential equations: 

ipo = -I22HQ - hzBfio, /3 = -I32H0 - l 33 B(3 (30) 

The homogeneous equilibrium stress is to — (c — Zi2&)"0o — li30o- We can see, 
that the coefficents of the above equations are nonnegative, therefore ipf = Q, 
/3f = 0, and tf = are asymptotic equilibrium values of the relaxation equations 

(20). 

4-2. Linearization 

In the following we assume small perturbation around the homogeneous 
equilibrium solutions, therefore the fields £ — (e, v, e, t, ip, /3) have the form 
£(t, x) — £o( T ) + S£( T , x). Then the linearization of the equations (fT9)) - (l23l) and 



(HHD-dlHJ) result in: 

d T 8e — d x Sv, (31) 

g d T Sv = djt, (32) 

g d T 5e = X ^^J dxxSe, (33) 

St = ESe + cStp + lud x 5v - gluibSif) + cSe - adxxSip) - Qh 3 B 5/3 (34) 

d T 5ip — hid x Sv — gl22(bSiJj + cSe — ad xx 5ip) — g^BSp, (35) 

d T Sf3 = l 31 d x 6v - gl 32 (b5ip + cSe - ad xx 6ip) - gl 2 3B5(3, (36) 



Here To is the temperature in the equilibrium and we have assumed that 
the temperature is the function of the internal energy and independent of the 
other state variables. Then (|33)l is independent of the other equations. Looking 
for conditions of stability we introduce exponential wave perturbations of the 
following form: 6£(t, x) = 8ac Tt+tkx . This exponential wave solution of the linear 
system of equations (|3TT) -(|36 ]) is damped if the real part of Y in the solutions is 
negative. The linearized balance of internal energy (1331) gives: 

gTSe = -A^pVfc. (37) 

Therefore Y = — A ^ d ^ a k 2 is negative because the Fourier heat conduction 

coefficient Xf/T 2 is positive and the specific heat d e T is positive to ensure non- 
negative entropy produstion and concave specific entropy function respectively. 

The remaining system of equations (|3Tj) . ([32"]) . (|34]) - (f36| is transformed to 
the following linear matrix form: 



/ 


r 


—ik 








° \ 


/6e\ 







QT 


—ik 








Sv 




-E + l 12 c 


—l\\ik 


1 


-c + l 12 (b + ak 2 ) 


-h 3 B 


St 




I22C 


—l 2 \ik 





Y + l 22 {b + ak 2 ) 


h 3 B 


Sip 


V 


I32C 


-hiik 





l 32 (b + ak 2 ) 


Y + l 33 Bj 


W 
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Exponentially growing wave solutions of ([55)1 emerge whenever V and k sat- 
isfy the dispersion relation obtained from the determinant of the matrix above 
where T has positive real parts. In order to simplify the calculations we assume 
that /13 = Z31 = and I12 — hi- In this case we preserve the viscosity In, the 
other dissipative terms in (|18p I22 and I12, and keeping also l 33 we have a term 
that is beyond the MES framework. 

The determinant is a fourth order polinomial in T: 

gT A + [g(bl 2 2 + Bl 33 ) + k 2 l n ] T 3 + 
[gbBA 23 + k 2 (E + bA 12 + 5/33*22)] T 2 + 
k 2 [l 22 6 + bBA + BEl 33 ] T + k 2 B5A 23 = 0. (39) 

Here A23 = '22*33 ~ *23*32 > 0, A12 = '22'n — *12*12 > 0, A = 'n'22'33 — 
'11*23^32 — '33*12*21 > are positive, because the symmetric part of lu is positive 
definite. We have introduced b = b+ak 2 and S = Eb — c 2 is nonnegative because 
of the convex free energy. Therefore the coefficients of the above polinomial are 
nonnegative. According to the Rout-Hurwitz criteria the real parts of a fourth 
order polinomial a^x 4 + a 3 x 3 + C12X 2 + a\x + ap = are negative, whenever the 
coefficients of the polinomial are positive and [16| 

Ci = 0302 — 0401 > 0, C2 — {a 3 a2 — a±ai)ai — aoa 2 > (40) 

In our case for (f3"9")l the first condition gives 

d= k A l 11 ( i i + Bl 11 l 33 + bA 12 ) + 

k 2 g(b 2 l 2 2A 12 + B 2 *ii* 33 *33 + * 22 (g 2 + 2bBl n l 33 )) + 

bBg(bl 2 2 + Bl 33 )A 23 > (41) 

One can see, that all terms are positive in the above expression. 
The second condition is more complicated: 

C 2 = 

k e l u [-B8l n A 23 + (fi+Bl n l 33 + feA 12 )(*22<5 + bBl u A 23 + Bl 33 (n-bl 2 12 ))] + 
fc 4 g [~2BSl n (bl22 + Bl 33 )A 23 + (6 2 *22A 23 + B 2 l n l 2 33 + l 22 B(B + 2bl n l 33 ))x 

x(-l 22 S + bBl n A 23 + Bl 33 {fi - bl 2 12 ))] + 
k 2 B 2 g 2 (bl22 + Bl 33 )A 23 (l 33 B 2 + b 2 A) (42) 

Some transformations may be usefull here. The coefficient of the k 6 term 
can be written as 

S 2 5 

-p(bl 2 2 + Bl 33 ) + - [b 2 (bl 22 + 2Bl 33 )A X2 + B 2 (Bl 33 + 2bl 22 ) + bB 2 l n l 2 33 ] + 

^B(B 2 + bB(A 12 + *ii* 33 ))(S 2 *33 + *> 2 A), (43) 
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and that is clearly positive. The coefficient of the k 4 terms follows as 

j(b 2 l 22 A 12 + B 2 l n l 2 33 + l 22 (B 2 + 2bBl n l 33 )) + 

- (Bl 33 + bl 22 ) [b 2 l 22 A 12 + l 22 B 2 + Bl u (Bl 2 33 + 2bl 23 l 32 ] (44) 

This is positive, too, except of the last parenthesis in the last term, where 
l 32 l 23 may be negative. It is reasonable to consider that term as a second order 
polinom of 6, which is positive, if the discriminant is negative. This condition 
is: 

Disc = l 22 {l 33 ) 2 A 12 - ln(l 23 l 32 ) 2 < 0. (45) 
This inequality can be rewritten in the following form: 

/3 3 A 23 A 12 + I 2 A + l n L 2 (2l 2 — L 2 ) > (46) 

Where we have introduced the I = (l 23 + 1 32 )/2 and L = (l 23 — l 32 )/2 for the 
symmetric and antisymmetric parts of the couplings of ij) and (3 evolution. A is 
the determinant of the symmetric part of lu 

A = lnl 22 l 33 - l 33 l 12 l 12 - lul 2 (47) 

The final condition (|46| cannot be violated if 2l 2 > L 2 , however the last term 
may dominate the previous ones for sufficiently large L. This can happen only 
if also or l 22 is not zero, otherwise the second term in (|4"4"]) is zero. Therefore 
dissipation is necessary for the violation of generic stability conditions. 



5. Discussion 

We have introduced a reasonably simple dissipative version of the microde- 
formation theory of Mindlin-Eringen-Suhubi. The generalization is based on 
thermodynamic principles extending the basic state space of continuum mechan- 
ics by dual weakly nonlocal internal variables. The most important distinctive 
property of the generalization is the nonclassical form of the entropy current 
density when the entropy density depends on the gradient of the state variables. 
There are two remarkable predictions of the new approach. The first one is the 
extra term in the thermodynamic force of heat conduction beyond the temper- 



ature gradient 12j. The second one is the modification of the nondissipative 
part of the stress. The second contribution is related the question of material 
time derivatives for the internal variables. A more detailed and more rigorous 
analysis of the thermodynamic approach to generalized mechanics based on the 
exploitation of the entropy inequality in case of second order weakly nonlocal 



constitutive state spaces in first Piola-Kirchhoff framework is given in [17 1. 

In a simple but reasonably representative example in one space dimension we 
have obtained the conditions of generic stability of the theory. We have shown 
that in some particular cases the thermodynamic conditions alone do not ensure 
generic statibility. 
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The stability of homogeneous equilibrium in case of small perturbations is 
a basic physical expectation in any dissipative continua, therefore we suggest 
generic stability as a benchmark for testing the viability of dissipative continuum 
theories both in classical and relativistic spacetimes. 
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